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ABSTRACT 

GRB X-ray light curves display rapid declines followed by a gradual steepening or plateau phase in 
> 30% of GRBs in the Swift sample. Treating the standard relativistic blastwave model in a uniform 
circumburst medium, it is shown that if GRBs accelerate ultra-high energy cosmic rays through a Fermi 
mechanism, then the hadronic component can be rapidly depleted by means of photopion processes on 
time scales ^ 10^ - 10'* s after the GRB explosion. While discharging the hadronic energy in the form 
of ultra-high energy cosmic ray neutrals and escaping cosmic-ray ions, the blast wave goes through 
a strongly radiative phase, causing the steep declines observed with Swift. Following the discharge, 
the blast wave recovers its adiabatic behavior, forming the observed plateaus or slow declines. These 
effects are illustrated by calculations of model bolometric light curves. The results show that steep 
X-ray declines and plateau features occur when GRB sources take place in rather dense media, with 
n > 10^ cm-3 out to > 10^"^ cm. 
Subject headings: gamma rays: bursts 
cosmic rays 



stars: winds, outflows — nonthermal radiation physics 



1. INTRODUCTION 

The Swift telescope is providing a new database of 
GRB light curves consisting of a BAT light curve in the 
15 - 150 keV range followed, after slewing within « 100 
s, by a detailed 0.3 - 10 keV XRT X -ray light curve 
and UVOT monitoring (jCehrels et al.l[2004f ). Extrapo- 
lating the BAT light curve to the XRT r ange gives ~ 1 
keV X -ray light curves since the trigger. lO'Brien et al.l 
()2006f ) present a catalog of the combined 0.3 - 10 keV 
light curves of 40 GRBs, of which 14 or so have mea- 
sured redshifts, w 30% display rapid X-ray declines, and 
an additional « 30% display features unlike simple blast 
wave model predictions. In some GRBs, the 0.3 - 10 keV 
flux can decrease by 4 or 5 orders of magnitude over a 
period of < 10^ seconds within several minutes after the 
GRB trigger (e.g., GRB 050915B, GRB 050422, GRB 
050819). About one- half the XRT sample shows X-ray 
flares or short timescale {At/t ^ 1) structure at > 10"^ 
s after the GRB trigger, and in some cases out to > 10^ 
s (e.g., GRB 050904 at z = 6.29 and GRB 050730 at 
z = 3.97). 

Following the rapid X-ray declines, a more gradual 
steepening commences in most XRT light curves, as in 
the cases of GRB 050219A or GRB 050607. In several 
GRBs, e.g., GRB050315 or GRB050822, a plateau phase 
with rising and decaying features within 0.1-1 day af- 
ter the GRB trigger are observed; this phase may also 
be found in other GRBs (GRB 050724 or GRB 050819), 
but could be contaminated or mimicked by a late-time 
X-ray flare. In a few GRBs (e.g., GRB 050401, GRB 
050717), the X-ray decline is gentle and monotonic, but 
in general the XRT light curves, whether displaying over- 
all convexity or concavity, reveal temporal structure and, 
oftentimes. X-ray flares (e.g., GRB 050712, GRB 050716, 
GRB 050726). A phenomenological model with two dis- 
tinct components can fit the rapid decays and hardenings 
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(IWillingal e et al.ll2006l ). but lacks a physical basis. An 
acceptable physical model for GRB afterglows must be 
able to explain this diversity of behaviors. 

A combined leptonic-hadronic GRB model is proposed 
in this paper as the cause of the rapid X-ray declines 
and plateaus discovered with Swift. The analysis pre- 
sented here fol lows the standard blastwave model (e.g., 
lMeszaro^ l2006\ We show that if GRBs accelerate cos- 
mic rays to ultra-high energies, then for certain classes 
of GRBs, the blastwave will become strongly radiative 
during the early afterglow, and consequently will exhibit 
rapid X-ray declines. This class of GRBs is defined by a 
range of blast wave and environmental parameters. One 
set of parameters that produces such declines has ini- 
tial Lorentz factors Fq ^ 100 - 300, apparent total en- 
ergy releases > 10^^ ergs (and absolute energy releases 
< 10^^ ergs), and surrounding medium density n > 10^ 
cm^"^, which is assumed to be proton dominated. The 
blast wave microphysical parameters eg, e_B must both be 
> 0.1. 

For these GRBs, Fermi processes in the blast wave are 
assumed to accelerate proton and ions, like they do elec- 
trons, to ultrarelativistic energies. By making reasonable 
approximations for the acceleration rate as a fraction of 
the Larmor rate, and particle escape through Bohm dif- 
fusion, we find that photohadronic losses and particle es- 
cape signiflcantly deplete the internal energy of the blast 
wave, causing the blast wave dynamics to be strongly 
affected. Photopion interactions by the ultrarelativistic 
protons and secondary neutrons with the internal syn- 
chrotron photons make a source of escaping neutrons, 
neutrinos and cascade 7 rays, in addition to a generally 
weaker proton synchrotron component. 

In Section 2, a standard blast-wave physics analysis 
for emissions from a GRB external shock in the prompt 
and early afterglow phase is presented, and timescales 
for the various processes are calculated for an adiabatic 
blast wave that decelerates by sweeping up material from 
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a uniform surrounding medium. Parameter sets that al- 
fow a large fraction of the internal energy to be radia- 
tively discharged through hadronic processes are graphi- 
cally examined in Section 3. In Section 4, the equations 
for blastwave evolution are solved in the case of inter- 
nal energy that is promptly radiated or is exponentially 
depleted with time. Synthetic bolometric light curves 
are calculated. Although these cannot be directly com- 
pared to the Swift X-ray light curves without taking into 
account spectral effects and more complicated radiation 
physics, the bolometric light curves exhibit many of the 
features observed with Swift. 

Discussion of multiwavelength and multichannel 7-ray, 
cosmic-ray, and neutrino predictions for this model is 
found in Section 5, including a comparison with other 
models for the X-ray declines. The study is summarized 
in Section 6. 

2. ANALYSIS 

The energy flux = dE/dAdt = vF,, = = 
L^/A'iTd\, where = 10^*^28 cm is the luminosity dis- 
tance, the source luminosity — Attx^cu'qT'^ , x is the 
distance of the blast wave from the explosion center, and 
r = V{x) is the blast wave Lorentz factor at x. The 
measured dimensionless photon energy e = hv/nieC^ is 
related to the emitted photon energy e' through the re- 
lation e' = (1 + z)e/2r, and the differential distance 
traveled by the blast wave during reception time dt is 
dx = 2T^cdt/{\ + z), where z is the source redshift, and 
primes denote quantities in the comoving fluid frame. 
Thus the comoving energy density 



where kkin ~ 1/3 gives the photon energy density that 
Toduces a rec eived spectrum with peak i/F^ flux = 
Dermeii 120041) . Relating x to the measured variability 
time gives the proper spectral density of the radiating 
fluid as a function of the Doppler factor. The GRB spec- 
trum is approximated by the broken power-law form 

where u = e/e^k = l^'pk^ fe^k = 10"*^ /-6 ergs cm^^ g-i 
is the peak vF^, flux at the vF^, peak photon energy epfc, 
a > and 6 < are the vF^ indices, and the Heaviside 
function H{y) — 1 for y > \ and H{y) = otherwise 
restricts the lower and upper branches of the spectrum 
to their respective ranges. 

2.1. Adiabatic Blast Wave 

Consider a blast wave with coasting Lorentz factor 
Taoo = ro/300 and apparent total isotropic energy re- 
lease Eo = 10^'*£'54 ergs, so that its absolute energy re- 
lease, due to collimation of the relativistic outflows, is 
wave sweeps through a uniform 
surrounding medium with proton density no = 100n2 
cm~^, it will slow down on the deceleration length scale 

* "m (3) 



= 2.6 X 10' 



(H eszaros &: ReesI fl993f) . A relativistic adiabatic blast 
wave decelerates according to the relation 

(4) 



r = ro/v/l + 2(a;/xd)3 



(jBottcher fc Dernieil |2000[) . from which can be derived 
the asymptotes 




and 



r < 1 



To |2-7/8r-3/8, r»l 



(5) 



(6) 



Here the dimensionless time r = t/td = t' /t'^, where the 
deceleration timescale is 

*rf = —F2— =9-6(1 + s, (7) 

i qC ll2^ 300 

and the inverse of the comoving deceleration timescale is 
t'-i = (^)-i ^ 3.5 X l0-4(!^^)i/3 s~i . (8) 

1 oC -t^54 

The available time in the comoving frame is 



T < 1 



i' — i' 

^ava ^d 



(21^8/5)^5/8 ^ o.872t5/8, t > 1 



, (9) 



2.2. Blast Wave Physics 

We tr eat the photopi on process in the fast cooling 
regime (jSari "eTallfTQQSh . The minimum Lorentz fac- 
tor 7mi„ = eeTprUpV /me, where Tp = {p - 2)/{p - 1) 
and 2 < p < 3. The emission detected by Swift is 
assumed to be predominately nonthermal synchrotron 
radiation. The mean magnetic fleld in the fluid frame 
is -B = bBcr ~ ^32TTmpc'^noeB^ — OAy/esno F G ^ 
370^eB-in2 r3oo(r/ro) G, where B^r = mlc^/eh = 
4.414 X lO^"^ G, and the minimum mean observed syn- 
chrotron photon energy from electrons with ^min is 

^ r67^„ _ o.85v^^rTTT^e2_i_^2^^rtoo(r/ro)4 



1 + 2: 1 4- z 

(10) 

The cooling Lorentz factor 7c — 3me(l + 
z)/lQmpa'TeBn{)T^ctdT, and the dimensionless cooling 
frequency (in units of rrieC^) is given by 

_ 2.3 X 10^(1 + z) _ 1.0 X io-«r3^o(r/ro)"^ 



(11) 

Comparing eqs. ([T0| and (fTT|) shows that we are in the 
strong cooling regime, ec < emin, when 



^-7.. in5/ T \2 4/3 i72/3t^8/3 

T < 7 X 10° (ee_ieB-iJ^5/2) n2 ^54 ^300 



(12) 



The photon energy at the peak of the vFjj synchrotron 
spectrum for the fast-cooling blast wave is given by 



0.85V^iZT7Tie2_i.F5/2rtoo / ^' 



r < 1 



l + z 



and the vF^, peak flux is given by 



(13) 



4 

(-CCTtJ/b) Iniin^'eil mm 

]x) . (14) 



47rd| '3 



In this expression, Ub — B'^/^tt is the magnetic- field 
energy density in the comoving frame. For the fast- 
cooling regime, N'^{^rnin) = ^e(a;)7c7mL' and Ne{x) = 
47rnoa;V3. Thus 



/,p,(10 ergs cm s 



2 ,-,-1^ 



2.7 



a: l/3p2/3Tn8/3 
(-e-l-r 5/2^2 -^54 ^ 300 



< 1 



(2r)-i, r>l 



(15) 



2.3. Photopion Losses 



The energy-loss rate due to photopion production on 
the GRB synchrotron radiation field is 



27;^' 



{^(i-y)[ 



2?/' 



a-l 



3 — a (3 — a)(a — 1) a— 1- 



^ max(l,y)'' 1 _ y'^mayi(l,yf ^ 
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after substituting eq. ^ into eq. p^ . with y = 
f'thrl'^l'p^'vk- H ere we use the approximation of 
lAtovan &: Permed ()2003[ ). where the product of the pho- 
topion cross section and inelasticity is CT^(eJ,) — a = 
70 /ib H{e'j. — ej^^), and the threshold dimensionless pho- 
ton energy for photopion production is ej^^ = 400 (i.e., 
ruec^e'th^ « 200 MeV). 
The asymptotes of eq. are 



,,6-1 



y > 1 

y<l, l<a<3 



}l^^^y-\ y«i,o<a<i 

(17) 

where y = 1 defines the Lorentz factor 7p of protons that 
interact primarily with internal synchrotron photons at 
the vF^ peak frequency e^j,. We call Ep ~ mpC^jp — 
nipC^T^p the peak cosmic-ray proton energy, as it is the 
characteristic energy of protons with Lorentz factor 7^ 
that would escape from the blast wave with Lorentz 7p 
factor as measured by a local observer. Hence the peak 
cosmic-ray proton energy is 



_ mpC^F^eJ,^,, _ 1.7 x 10i6(r/300)2 



and 



^^^077 (7p 



2r2mec2e' (l-6)(3-6) 



1.1 X 10-6fcfc^„d2g/_ 



pk 
28 /-6 



(l-6)(3-6)(l + z)a;2g(r/300)epfe 
where a; = lO^^xig cm. 



(19) 



2.4. Rates and Limits for Ultrarelativistic Protons 

2.4.1. Adiabatic Loss Rate 

Adiabatic expansion losses operate on the same 
timescale as the available time, so for comparison of the 
adiabatic loss rate to other rates, we write t'~^l = 
using eq. 

2.4.2. Photopion Loss Rate 

Substituting eqs. dS]), ([6]), and (fT3| into eq. p?|l gives 
the comoving photopion energy-loss rate 

5.2 X IQ-'^kkrnnl^^ s-i [ 1' 

(1 _ 5)(3 - 6)ri„oe^/!iee-l \ 223/8^3/8^ ^ » 1 



at the characteristic energy 
4 X 10^^ 



(20) 



< 1 



£^p(eV) 



VeB-l"2ee_i^5/2r3oo [ 3.4r3/4^ T > 1 



. (21) 

of an escaping cosmic ray, as measured by a stationary 
observer in the local source frame. 

2.4.3. Acceleration Rate 

Because a significant energy gain by a particle can take 
place through Fermi acceleration mechanisms on times 
not shorter than the Larmor time t'^^ = mc^'yp/eBc = 
[mc/ eB)(-f /T) (jRachen fc MeszaroslfTOOSl) . the accelera- 
tion rate in the proper frame can be written as race — 
^'acc — Cacct'{^^ , with the acceleration parameter Qacc ^5 1- 
Hence the acceleration rate at the peak cosmic-ray pro- 
ton energy Ep is given by 



^acci.Bp') — 



25Cocc n2e_B_lF3ooee_iJ^5/2 s ^ 



1 



11.4r3/2 ' 



T < 1 



r > 1 



(22) 

The acceleration rate for 102'^i?2o eV cosmic ray protons 
is 

racc{E2o) = 

1.1 X 10-=^(Cacc/0.1)VeB-i"2 F§oo _i 



E20 



(1 -f 27/8r3/8)2 



(23) 

Here we consider a standard acceleration parameter 
(acc = 0.1. Values of (acc — 1 require unreasonably effi- 
cient particle acceleration, with particles gaining a large 
fraction of their energy in a single Larmor timescale. 
If Cacc ^ 0.1, then GRBs would not accelerate cosmic 
rays sufficiently rapidly to make ultra-high energy cos- 
mic rays. 

2.4.4. Escape Rate 

The mean escape rate using the Bohm diffusion ap- 
proximation is given by ig^,^ — {x)"^ /2k-q, where kb — 
c^t'^/i is the diffusion coefficient. For particle acceler- 
ation in GRB blast waves, the characteristic dimension 
{x) is the shch width A' = /a^/F, and /a = 1/12 for 
the width of the shocked ffuid shell swept up from the 
circumburst medium by an adiabatic relativistic blast 
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waves fe.g.. iPanaitescu &: Meszaroill999|) . Thus the es- 
cape rate is 



''esc 



1 



(24) 



where C,esc is a parameter that allows particle escape on 
timescales shorter {Qsc > 1) or longer (Cesc < 1) than 
the escape timescale set by Bohm diffusion. If the GRB 
blast-wave shell entrains a randomly oriented, tangled 
magnetic field, then depending on the coherence length 
of the disordered magnetic field, the particles could dif- 
fuse more rapidly than given by Bohm diffusion, so that 
Cesc ^ 1- In contrast, if the GRB blast wave is assumed 
to entrain an ordered field, for example, a toroidal geom- 
etry in a jetted fireball, then escape could be impeded 
compared to the Bohm timescale, so that (esc 1. Here 
we consider the Bohm limit for the rate at which parti- 
cles escape, keeping in mind that the actual escape rate 
could be quite different. 
The escape rate for protons with characteristic energy 
IB), is 



Ep, eq. 



5 X lO-^Ce^c 



< 1 



1/3 2 T- T^2/3 I -I 

es-mj' ee-i-^5/2r3()o t 2At^/^, t > 1 



(25) 



The escape rate for 10^°i?2o eV cosmic ray protons is 

^e7c (-^20) = 



1 V in-3F r,l/'5p4/3> 1 f r<l 



/7 p2/3 



(2t)-i/2, r»l 



(26) 



Ef{10^° eV) = 2.4 e 



1/2 ^l/6pl/3pl/3 
B-l'h ^ 300^54 



25/8 , 



2.4.5. Size Scale Limitation 

We also have the iHillasI (|1984f l condition that the Lar- 
mor radius be smaller than the characteristic size scale 
of the system, which is the shell width A'. Requiring 
tl — {mc'^/eB){j/r) < A' = f/\x/T implies a limit to 
maximum proton energy, given by 

T > 1 ' 
(27) 

using the asymptotes, eqs. ([5]) and ([6]). Note the slow 
late-time decline oc t' '^/^ when r > 1 ()Vietril 119981 : 
iBottcher fc Dermeif[T998[ ). Thus we see that standard 
parameter values allow Fermi acceleration of protons to 
ultra- high energies in GRB blast waves when Cacc ^0.1 
and C,esc — 1, making GRBs a viable candidate for 
UHECR production. 

2.4.6. Proton Synchrotron Energy Loss Rate 

The inverse of the synchrotron energy-loss timescale 
for an escaping proton with Lorentz factor 7^ = T^'^ is 



I p,syn — ''p,syn 



ip,syn I 



16 ,me.2 

3 "-nip' 

^ 9.4 X 10-%ei3-ir3ooS20 (^) s^^ . (28) 

The mean proton synchrotron photon energy, in units 
of rrieC^, from protons with energies 10^°£'2o eV (£'20 — 



TABLE 1 

Parameter Sets Used in Analysis 





Stnd. Sot 


Set 1 


Set 2 




z 


1 


1 


1 




1 


1 


0.5 




1 


1 


10 




1 


10 


10 


ee-l 


1 


3 


1 




1 


3 


3 


Cacc 


0.1 


0.1 


0.1 


Cesc 


1 


1 


1 




19.3 


8.9 


122 


xa(W^^ cm) 


2.6 


1.2 


4.1 



^ Apparent isotropic energy release in 
units of 10^** ergs.'' Circumburst medium 
density, assumed uniform and dominated 
by H, in units of 100 cm"''. 



Standard Parameters 



(a) Calculated at E 

p 

'*^. E (GeV) 




10^ 

10' 

m 

in° 

10-' f 
10"^ oj 
10-3 a. 

1 X 
CD 

10= " 
10 = 
10' 



Fig. 1. — (a) Characteristic rates and energies calculated at 
Ep, (b) characteristic rates calculated at Ep = 10^'^ eV, using the 
Standard Parameter set. 

7ii = 7p/10^^) as measured in the stationary frame, is 
given by 



3r B m„ 



'^p^syn 



— ^4.5x10^ 



1 + Z Brr m. 



eB-l"2£^20 

l + z 



(29) 

independent of time in the relativistic deceleration 
phase — provided of course that e_B is time- independent 
and rt2 — const. 

3. RESULTS 

For initial parameter estimation, we adopt the Stan- 
dard Parameter Set given in Table 1, with z = l,r3oo = 
l,i?54 = l,n2 — l,ee_i — 1, and €3-1 — 1- This 
set is motivated by values that reproduce typical peak 
fluxes and durations fo r BATSE GRBs at z sa 1 
(jChiang fc Dermeil I1999D . except that here cb = 0.1 
rather than 65 = 10~^. 
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10-= 
10' 
10' 



Parameter Set 1 



Parameter Set 2 



. ^ (a) Calculated at 




E (10'^eV) 1 






E (GeV) 
^ r 

JS£ 






r" 1 

p,syn 


10 100 


1000 


t(s) 





10" r 



10" m 

10"' ? 



10" ? 




Fig. 2. — As in Fig. 1, but calculated for Parameter Set 1, with 
the energy of electrons, protons, and magnetic field in approximate 
equipartition. 

Fig. 1(a) shows the Standard Parameter rates of ac- 
celeration, photopion losses, proton synchrotron losses 
and escape for cosmic rays with energy Ep, the i^Fi, 
peak photon energy epk, the cosmic-ray peak energy 
Ep, and the mean proton synchrotron photon energy 
Ep^syn — 'meC^£p,syn radiated by protons with energy 
Ep. Here and throughout we use an acceleration fac- 
tor Cacc = 0-1 5 escape factor Cesc = Ij f^nd kinematic 
factor kkin — 1. 

From the top panel in Fig. 1, one sees that the accelera- 
tion rate exceeds the inverse of available time throughout 
the_early afterglow phase, so cosmic rays with energies 
Ep are in principle easily accelerated through Fermi 
processes to energies > Ep. Only at several hours into 
the afterglow do photopion losses limit acceleration to 
Ep < Ep, which by then is > 10^^ eV. At these late 
times, proton synchrotron emissions make a ^1% contri- 
bution to the total loss rate. The diffusive escape rate 
of protons can appear as a >1% effect on the total rate, 
but is generally insignificant in the early afterglow. The 
external shock emission from this GRB is brightest « 20 
s after first being detected, though shell collisions could 
make brighter features during the prompt and afterglow 
phases (see Discussion). 

The acceleration, escape, and loss rates for a cosmic 
ray proton with Ep = lO^*' eV are plotted in Fig. lb. As 
can be seen, there isn't enough time to accelerate cos- 
mic rays to > 10^" eV energies for these parameters, so 
there cannot be significant > 10^° eV (supcr-GZK) cos- 
mic ray production or photopion losses from such GRBs 
(unless (acc ^ 0.1). Protons would also escape before a 
significant fraction could be accelerated to such energies. 

A set of parameters that overcomes these limitations is 
easily found. Consider Parameter Set 1 in Table 1, with 
z = l,ro = 300, i?54 = 1,7^2 = 10, ee_i = 3, es-i = 3, 
giving the rates, fluxes, and energies shown in Fig. 2. For 



..^ (a) Calculated at 



10" 

10' 

10° m 

10"' ? 
tg 

10"^ to' 

10"' s 




10 100 1000 

t(s) 

Fig. 3. — As in Fig. 1, but calculated for Parameter Set 2. 



10^° eV cosmic rays shown in the lower panel, an interest- 
ing conjunction occurs when r^j^ « l/Cua ^ i^acc, which 
happens here at « 300 s. Protons accelerated to ~ 10^° 
eV energies are converted, ~ 1/2 of the time^, to ultrarel- 
ativistic neutrons that escape from the blast wave to form 
one c omponent of a neutral beam ijAtovan &: DermeH 
l2003f ). in addition to neutrinos and 7 rays. 

The GRB formed in Parameter Set 1, Fig. 2, has a 
rather high Epk ^ 8 MeV — which may be irrelevant in 
an internal/external scenario — but values of Fq < 300 
will lower epk during the prompt phase and lengthen the 
prompt phase duration. Fig. 3 shows the results for Pa- 
rameter Set 2 with z — 1, Fq = 150, Eq — 10^^ ergs, n2 = 
10, Ee-i — l,e_B_i = 3 . This GRB peaks > 50 s after 
the trigger, has a lower epk, and reaches a slightly lower 
/e^^ peak flux than in Fig. 2. 

Parameter Sets 1 and 2 model fast-cooling GRBs with 
eg = 0.3 and eg = 0.1, respectively, that exhibit a radia- 
tive photopion phase. By letting ee= 30% for Set 1, it is 
understood that a large fraction of the swept-up power 
is found in nonthermal electrons rather than in baryons 
or fields. A large body of parameter values clustering 
around the Parameter Set 1 values predict strong pho- 
topion losses in the early afterglow phase, especially if 
Cacc ^ 0.1 is allowed. If eg ^ 0.1 is assumed, as in Pa- 
rameter Set 2, then agreement with GRB energetics for 
GRBs with apparent bolometric j-ray energy > 10^"* ergs, 
which represents a signific ant fraction of pre-Swift GRBs 
(|Friedman fc Blooml l2005l ). means that a large fraction 
of GRBs also have apparent total energy E^^ > 10, due 
to the unknown efficiency of converting the total energy 
into 7 rays. Even for small ee models, the necessary large 
values of total apparent total energy can result in strong 

^ The charge-changing fraction is not 1/3, as expected for A 
resonance excitation and decay, due to the inclusion of direct pion 
channels above threshold, and multi-pion production at energies 
far above threshold. 
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photopion losses. Thus if GRBs accelerate UHECRs, 
then photohadronic losses efhciently deplete GRB inter- 
nal energy and affect the blastwave dynamics over a wide 
range of GRB parameters. 

4. BLAST WAVE EVOLUTION WITH RADIATIVE 
DISCHARGE 

The previous section showed that with reasonable pa- 
rameter values, photopion production and escape can 
rapidly deplete internal blastwave energy in the early af- 
terglow when hadronic processes become important. We 
now calculate model light curves resulting from the radia- 
tive discharge of internal energy through photohadronic 
processes. The dynamics of the blast wave is numeri- 
cally solving using the equation of relativistic blastwave 
evolution, including adiabatic losses, given by 



dP/d. _TP^ + {^^^ 



r2 - 1 Mo + m{x) + U' 

(jPermer fc Humill200lD . Here the blast wave momentum 

P = vT^ — 1, Mq — Eq/To is the initial baryon loading, 
that is, the baryonic mass mixed into the explosion, and 
m{x) is the swept-up mass (all masses are now in energy 
units). The internal energy, excluding rest mass energy, 
is given by 

/"OO 

[/' ^ U'{x) = mp / dp'ijL - l)N'ip'; x) = 



dxi xf no{xi) (7 - 1) , 



(31) 



where p' = — 1 is the proton's dimcnsionless mo- 

mentum in the comoving frame, 7^ its Lorentz factor, 
and 

p = ^,) = = Pi^^) (-) [^] (32) 

solves the equation for momentum evolution of a non- 
thermal proton in a blast wave with thickness = 
fAx/ T ()Panaitescu fc Meszarod Il999t [D ermer fc Humil 
I2OOID . Here only adiabatic losses are considered for the 
protons, and /a is assumed to remain constant with x 
and take the value 1/12. The differential swept-up mass 



dm{x) 
dx 



AirnipX^ no{x) , 



(33) 



where no{x) is the circumburst medium density, assumed 
to be radially symmetric about the GRB source. The 
internal energy U' changes due to volume expansion and 
adiabatic losses according to the relation 



dULu 
dx 



^(1 



J dx, xf no{x^){'^) . 

We consider two case scenarios to simulate the change 
in internal energy due to a photohadronic discharge: 

1. A fraction ^ of the internal energy is instanta- 
neously radiated away at x = Xq. Following the 
discharge, the blast wave is again assumed to sweep 
up circumburst material and evolve adiabatically. 



In this case, U'{x) and dU'^^^/dx are given by eqs. 
([3T|) and (l34|) . respectively, a,t x < xq. At x > xq, 



(35) 



^ -(1-0/ dx, x^ noix,) i-f - l) + 



Airnip jQ 

dxi xf no{xi) (7 - 1) 

Jxq 

and 

dx X 3 dlnx 

rxo p2 rx -2 

[(1-C) / dx, xf no{xi){—)+ dxi xf no{xi){—)]. 

1 J XQ 7 

(36) 







2. The internal energy is radiated on an exponential 
dissipation timescale t'^,,^ , followed by a recovery of 
the blast wave to adiabatic behavior at 2; > xq. 
Thus 



U'{x) 



Attitl 



dxi xfno{xi){'^ - l)exp[~t'{xi)/ti,^t] + 



p Jo 

H{x ~ Xq) I dx, xj no{xi) (7 - 1) , (37) 

•'Xo 

and 

dKdi ^ ^TTTOp _ 1 dlnV ^ 
dx X 3 d\nx 

[ [ dx, x\ nQ[x,){^) exp[-i'(a;,)/t^„t] 
+ H{x-xo) / dx, xf no{xi){—)] , (38) 

Jxo 7 

for all a; > 0. 

We consider first Case 1 with Parameter Set 1,^ with 
a sudden loss of radiative energy taking place at xq — 
5 X 10^^ cm, corresponding to an observer time of w 
650 s after the start of the GRB — essentially the same 
timescale when photopion losses and escape become im- 
portant (Fig. 2b). In Fig. 4(a), the numerical solutions 
to the equations for the evolution of P are plotted for 
^ = 0,0.5,0.9,0.99. Also shown are some approxima- 
tion s to the dynamics of the blast wave when ^ = 
(seeE ermer fc Humil [2OOI : the curve labeled "momen- 
tum conservation solution" has the term dU'^^^Jdx set 
equal to zero; the curve labeled "analytic blast wave" 
is eq. the curve labeled "momentum conservation 
solution" s hows the analytic form for a fully radiative 
blast wave (|Blandford fc McKeel[T976 l: lChiang fc Permed 
119991 ). The inset gives P as a function of observer time 
using the numerical solution to eqs. ((30| . p5|) . and ((38| . 
for different values of C,. 

Fig. 4b shows the terms PT^dm/dx) and ~{P + 
P~^){dU' /dx)adi in the numerator of eq. ([501) fo'' P^" 
rameter Set 1 with ^ = 0,0.5,0.9, and 0.99, and Fig. 

^ The microphysical parameters are irrelevant in this treatment 
of blastwave dynamics, though this may not be true in general, 
particularly for radiative GRBs with es-i > 1 in the fast cooling 
regime. 
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x(cm) ^jc^j 

Fig. 4. — (a) Solid curves show the evolution of a GRB blast wave 
as a function of x for different fractional discharges 5 = 0, 0.5, 0.9, 
and 0.99, for Parameter Set 1. Also shown are approximations to 
the case with ^ = (see text). Inset shows blastwave evolution as 
a function of observer time, (b) Variation of values of the terms in 
the numerator of eq. lijOt . as labeled, (c) Variation of values of the 
terms in the denominator of eq. II30I I. as labeled. 

4c shows the terms Mq, m[x) and U'{x) in the denomi- 
nator of eq. ([SO]) . A reasonable approximation to blast- 
wave dynamics can be obtained, when the blast wave is 
relativistic, by neglecting the {dU' /dx)adi term because 
this term is a small, constant fraction of the sweep-up 
term PT(dm/dx). When the blast wave becomes non- 
relativistic, at a; = 4 x 10^^ cm for these parameters, 
the {dU' /dx)adi term must be retained. In Fig. 4c, the 
range of x where the Mq, U'{x) amd m{x) terms domi- 
nate the value of the denominator of eq. ([50)1 define the 
coasting, relativistic self-similar, and non-relativistic self- 
similar regimes, respectively. In the Case 1 scenario, after 
the instantaneous discharge at xq, the blast wave rapidly 
decelerates, but then recovers its adiabatic behavior as 
it sweeps up additional material from the circumburst 
medium. 

An instantaneous discharge may be oversimplified, 
however, because a discharge of any sort must take place 
over a finite time exceeding at least Ar'/c, the light travel 
timescale across the blastwave width. For this reason we 
now consider the Case 2 Scenario, again with Parameter 
Set 1. Here the blastwave is assumed to recover its adi- 
abatic behavior at = 5 x 10^^ cm after suffering an 
exponential discharge with comoving (exponential-)loss 
timcscales t'^^^ = 10, 20, 40, and 80 ks. Although this is 
sufficient to solve the blastwave dynamics, at least in the 
limit ee_i < 1, it still is not clear what is the shape of 
the light curve. 

We can derive some idealized light curves by weight- 
ing internal energies and rates by the beaming factor 
that boosts energy and rate in a spherical blastwave. 
The following weightings are employed, keeping in mind 
that this sort of approach neglects particle cooling and 
spectral effects. The weightings considered are 

1. r^C/', so that the comoving power is assumed to be 
proportional to the total internal energy, as might 
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Fig. 5. — Idealized light curves for the Case 2, Parameter Set 1 
Scenario with = 5 X 10^® cm and t'^^^ = 10, 20, 40, and 80 ks. 

correspond to a single species emitter that radiates 
in the adiabatic limit; 

2. r^L', where the comoving power V is set equal 
to the swept- up power A-Kx^nQnipC^ (3{T'^ ~ T), so 
this case would represent the bolometric leptonic 
luminosity of a GRB radiating in the strong cooling 
regime; 

3. r^C/'^, so that the comoving power is proportional 
to the square of the internal energy, as might hold 
in a scenario involving wave-particle coupling; and 

4. V'^U''^ /{x/xdY , where the comoving power is as- 
sumed proportional to the product of the total in- 
ternal energy and energy density, for example, of 
the magnetic field. 

Fig. 5 shows some numerical solutions for the Case 2, 
Parameter Set 1 Scenario with the different weightings 
just described. The blast wave is assumed to recover 
its adiabatic character at x^/xd = 4.2. The notation 
here is that U' = 10*C/^ ergs and L' — lO^L^ ergs s~^. 
Overlooking the shape of the light curves at early times — 
a point to which we return soon — these synthetic light 
curves resemble the Swift BAT/XRT X-ray light curves 
of GRBs. On top of these generic light curve shapes 
are the X-ray flares made by some as yet unspecified 
mechanism. 

The various weightings reflect different underlying as- 
sumptions of the physical model and yield a variety of 
temporal behaviors that could explain the range of X- 
ray light curves observed with Swift. The sensitivity of 
the model light curves to the location xq where the blast 
wave begins to evolve adiabatically and the comoving dis- 
charge timescales is illustrated in Fig. 6, showing model 
light curves with a range of X-ray declines and plateau- 
like features. 

5. DISCUSSION 

The central idea proposed here is that the rapid X- 
ray declines and plateaus discovered in ~ keV light 
curves of GRBs by the Swift team (jBarthelmv et al.) 
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Fig. 6. — Idealized light curves for the Case 2, Parameter Set 1 
Scenario with (a) zq = 3 X 10^^ cm and t'^^^ = 2, 4, 8, and 16 ks, 
and (b) xo = 6 x 10^^ cm and t'^^^ = 20, 30, 40, and 80 ks. 



l2005HTagliaferri et al.l [2005 HO'Brieii et al.ll2006D are sig- 
natures of UHECR acceleration in GRBs. The solu- 
tion to the problem of UHECR origin may involve a 
variety of source classes, but must involve at least one 
source class, and GRBs offer a very attractive solution: 
They are foun d outside our Gal a xy (and possibly i n 
our Galaxy; seelWick et all (|2004D : iMelott et all (|2004D : 
iDermer fc Holmes! (|2005lV r which must be the case to 
make UHECRs with Larmor radii larger than can be con- 
tained in the Milky Way; they are energetic, with each 
explosion releasing as much as ~ 10^^ ergs kine tic energy 
([Friedman fc Bloom|[2005l iLe fc Dermeij|2007l ). much of 
it in very clean, highly relativistic bulk particle or field- 
dominated outflows; and they are impulsive, at least on 
recurrence times > day. This last is important by lead- 
ing to a prediction of no clustering or angular corre- 
lations on the sky of UHECRs or iginating from GRBs 
(jWaxman fc Miralda-Escuddll996h . 

Simple (pre-BeppoSAX) BATSE estimates for the re- 
quired local emissivity in super-GZK (> 10^° eV) UHE- 
CRs compared with the emissivity in X-rays and 7 rays 
from GRBs, assumed to be at a mean redshift (z) sa 1, 
agree to within ~ 1 order of magnitude, suggesting a 
connection be tween the two (|Waxmanlll995l : IVictri 1995; 
lDermeil[200l ). If GRBs are the progenitors of UHECRs, 
then the superposed cosmic-ray intensity spectrum from 
all the UHECR GRB sources over cosmic time display 
features of propagation. Most pronounced are the GZK 
feature and the ankle fea ture, the latter of which is likely 
due to pai r production ("Berezinskii fc Grigo r'eval 119881 : 
I Wick et al . 2004; Bcrczinsky ot al. 2005). li Auger mea- 
sures a spectrum similar to the HighRes spectrum at en- 
ergies > 10^^'^ eV, then these spectral features of cosmic 
rays, accumulated during transport from the source to 
the detector, are reasonably well-explained within a sce- 
nario where UHECRs originate from cosmic rays, with 
clearcut predictions for GZK neutrino fluxes. 

Here we extend the scenario that UHECRs originate 
from GRBs by entertaining the hypothesis that the X- 
ray light curves of GRBs show a spectral feature associ- 
ated with UHECR acceleration. The synthetic bolomet- 
ric light curves in Figs. 5 and 6 are idealized, but contain 
the basic features of more realistic calculations that de- 
pend on a wide range of underlying assumptions about 



particle acceleration and the fraction of energy going into 
magnetic fields and waves. Depending on the radiation 
model, the internal and external shock contribution, and 
the setting of the zero of time, a wide range of X-ray 
decays can be understood within this picture. 

An important feature of this scenario is that a large 
fraction of the swept-up internal energy resides in pro- 
tons accelerated to ultra-high energies. Photopion 
processes effectively discharge the particle energy in 
the form of a ne utral beam (jAtovan fc Dermeil l2003l : 
iDermer fc Atova"n|[2004.) . The decay of the internal par- 
ticle energy density weakens the magnetic field through 
feedback of the particle energy into field energy, halt- 
ing further acceleration while hastening escape of the re- 
maining UHECRs found in the blast wave, until a dom- 
inant fraction, > 90% (depending on the initial amount 
of energy contained in the nonthermal protons), is lost 
from the blast wave as a neutral beam. 

Developing detailed physical models to simulate lep- 
tonic and hadronic acceleration and loss is beyond the 
scope of a single paper, but the various physics issues 
that enter into a combined leptonic/hadronic GRB model 
can be described in somewhat more detail. We break the 
discussion into a consideration of 

1. the acceleration mechanism; 

2. light curves in an internal/external shock scenario; 

3. light curves in an external shock scenario; 

4. explanation for the rapid X-ray declines; 

5. explanation for the plateau phase; and 

6. predictions of this scenario. 

5.1. Acceleration 

Acceleration to super-GZK energies is possible 
via stochastic proce sses in the blast-wave shell 
(IDermer fc Humi l l200lD and through acceleration by the 



internal shocks in a colliding shell model (e.g., Waxma3 
Il995l) . Acceleration to such energies through shock Fermi 
processes is not possible through relativistic external 
shocks forme d in a surrounding medium w ith magnetic 
field ~ ^G (jGallant and Achterberj 11999 1 but might 
be possible if the GRB takes place in a highly magne- 
tized environment, as might be expected if the surround- 
ings a.re formed by the stellar w inds of high-mass stars 
(e.g., IVolk and BiermannI Il988l noting limitations im- 
posed by total energy contained in the magnetized wind) , 
or if the upstream field is ampli fied by strea ming cosmic 
rays, for example, by thelBell fc iTuceld ()2001[ ) mechanism 
(|0'C. Drurv et al.l[200l . 

Here, acceleration of swept-up protons to ultra-high 
energies is assumed to take place by second-order gy- 
roresonant processes in the blast wave shell. This is 
reflected in the derivation of the acceleration rate, eq. 
(|23|) . which contains only shocked fluid quantities. A 
turbulent, stochastic Type 2 Fermi mechanism in the 
blast-wave shell formed by a relativistic (internal or 
exter nal) shock can make a highly efhcient accelera- 
tor ( Dorm er fc Humil l200lV whi ch is a sort of turbu- 
lent b oiler discussed decades ago (jOchelkov fc Prilutskiil 
Il975l and references therein), though here found in 
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the shocked fluid sheU of a relativistic blast wave. 
The basic coupling involves gyroresonant acceleration 
of ions and electrons interacting with MHD wave 
turbulence (in the context of Solar flares, see , e.g. , 
gteinackcr and Miller 1992"; iMiller a nd Ramatv 1989; 
Miller and Roberts _1_995) , the form of which is model- 
dependent and may furthermore involve anisotrop ic cou- 
pling depending on wave type (|Goldreichl 1200 1[ ) . But 
for isotropic power-law wave turbulence, acceleration 
through the stochastic Fermi mechanism makes hard 
number spectra (^(7) oc 7"**), with s « 1 and most 
of the energ y content consequently i n the highest parti- 
cle energies (lSchlickeisei] ll984 [l989l : iDermer et aLlfl99fil : 
iDermer fc Humil " 20011 ) . Combined first - and s e cond- 
order processes (e.g., I Schhckciser TgP; iKriillsl [l99l 
lOstrowski fc Schlickeiser 1993) suggest that steeper 
spectra, with 1.5 < s < 2, would be formed in this sort 
of scenario. But if s < 2, most of the energy resides in 
the highest energy particles, which can be accelerated to 
energies > 10^0 eV ()Dermerll2006l ). 

The gyroresonant wave-particle interactions acceler- 
ate particles dynamically, so no steady state is reached 
(analytic solutions to time-dependent particle evolu- 
tion through seco nd-orde r processes are given by, e.g . , 
iPark fc PetrosianI 119951 iBecker. Le. fc Dermeil 120061 ). 
and energy flows between waves, particles, and fields un- 
til it is discharged from the system. The bulk kinetic en- 
ergy of the blast wave, initially dissipated in the shell as 
highly nonthermal internal particle kinetic energy as well 
as field and wave energy, is transformed into a compo- 
nent of ultra-relativistic protons carrying a large fraction 
of the total energy. This internal energy, including the 
wave energy that continues to be fed into the particles, 
is discharged when photopion losses become sufficiently 
great. Growth of the target photon field, for example, 
from proton synchrotron and photohadronic secondaries 
followed by electromagnetic cascading, can lead to an 
photohadronic or proton -synchrotron loop ins t abilit y of 
the type considered by iKirk fc Mastichiadii (|1992D to 
extract increasingly more energy of the UHECRs until 
the internal energy content is effective discharged. Ob- 
viously, more model studies are required to quantify the 
total fraction of internal energy that can be discharged 
in this way. 

Under the given circumstances, ~ 10^® eV photopion 
neutrinos are created from the decay of photomeson 
secondaries, ultra-high energy neutrons escape if they 
avoid further photopion interactions, and an electromag- 
netic channel consisting of 7 rays and e"*" and e~ ini- 
tiates a cas cade, ultimately lea ding to a flux of escap- 
mg 7 rays (|Atovan fc Permed [20031. The cosmic-ray 
neutron discharge from GRB blast waves, as considered 
here, produces a component of the UHECRs, causing the 
ionic composition of UHECRs to be proton-dominated. 
The decaying UHE neutrons and, indeed, photoprocesses 
by the highest energy p rotons and ions make an UHE 
synchrotron decay halo (DermeJ 120021 : lloka et all l2004( l 
around host galaxies. The mean injection spectrum per 
GRB by this process has not been calculated from first 
principles, but could range from a fiat 7^71(7) over a nar- 
ro w energy range ^ 10^^ - 10^^ eV, as in the model 
of IWaxman fc Bahc alf (1999), to a power-law injec tion 
7^n(7) oc 7~°-^, as in the model bv lWick et al.l (|2004 ). In 
principle, this process could also operate in radio/7-loud 



AGNs, and produce s pectra as steep as "Pjiij) oc 7^"- 
found in the models of iBerezinskv et aT] (|2006f l. 

The acceleration mechanism in a colliding shell sce- 
nario is much different than stochastic acceleration in 
the fiuid shell, but the photohadronic discharge mecha- 
nisom could certainly operate in such a system, though 
the details still need to be worked out. Incidentally, the 
potential importance of second-order pr ocesses in GRB 
blast waves may negate the concerns of i Ghisellini et all 
(2000) that the GRB must display a strong-cooling spec- 
trum. Second-order acceleration competing with lep- 
tonic radiative and adiabatic energy losses can form a 
low-energy pileup in the electron distribution, so that 
no cooling tail appears. Indeed, the shape of broadb and 
time-averaged spectra of GRBs ( Schaefer et al.|[r998D . in- 
terpreted as nonthermal leptonic synchrotron radiation, 
seems to require an abrupt low-energy cutoff and may 
be consistent with an electron pileup at some lower en- 
ergy due to combined first- and second-order processes 
with cooling. This question will take greater urgency 
when GLAST provides combined GLAST GBM, GLAST 
LAT, and Swift BAT and XRT GRB light curves, per- 
mitting detailed spectral analysis of bright GRBs from 
the prompt to afterglow phases. This will reveal the tran- 
sition from the internal shocks (in the internal/external 
shock scenario) to external shocks, as now considered. 

5.2. Internal/External Model 

In the intern al/external scenario (|Pira^[l999L [2OO5I : 
iMeszaro^ |2006[ ) . the transition from the prompt phase 
formed by colliding shells of jet plasma to the afterglow 
phase made by external shocks takes place no later than 
the time of BeppoSAX reorientation, ~ 8 hrs, and likely 
much sooner. In this picture, the internal shocks make 
a distinct radiation signature from the external shock 
emission received later. The zero of time for the ex- 
ternal shock emissions from the GRB can, in the inter- 
nal/external scenario, be set to various times, but is de- 
termined observationally when the GRB detector is trig- 
gered, e.g., by precursor emission. In the case of a GRB 
exploding in a uniform circumburst medium, the latest 
the zero of time can be set would be about when the ex- 
ternal shock signature reaches its maximum brightness 
and triggers the GRB detector, which is generally near 
the dece leration time td- T he freedom to set the zero 
of tim e (jLazzati fc Begelma n 2006; Kobayashi fc ZhanJ 
l2007f ) when modeling data with the external shock com- 
ponent considered here, as shown in Fig. 7, adds an addi- 
tional distortion and increased variety to the model light 
curves, similar to various behaviors of the X-ray light 
curves found with Swift. 

An alternate view of the internal/external model is 
to have all the prompt emissions up to the time of the 
rapid X-ray declines made by internal shock processes 
(see §5.5). Here it is surprising that the Swift BAT 
light curves are rather smooth and apparently not a su- 
perposition of generic kinematic cur vature light pulses 
(jYamazaki et al.ll2006t rDermeill2004l ). although the last 
pulse in some GRBs, for example. Sw ift XRT data of 
GRB 061202 ~ 120 s after burst trigger ()GCN Report 19l 
l2006( l. do appear possibly to represent the last major ac- 
cretion events from a hypothetical coUapsar torus. In 
either variant of the internal/external scenario, a com- 
plete model GRB light curve would combine the model 
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Fig. 7. — Idealized light curves for the Case 2, Parameter Set 1 
Scenario with xq = 5 X 10^^ cm and i^^^ = 10 ks, with the time 
reset to -3 s and + 10 s for the thin and thick curves, respectively. 

light curves calculated in an external shock treatment 
with a second component consisting of 7-ray pulses and 
X-ray flares from colliding shells made by central engine 
activity. 

Placing these calculations in the context of the inter- 
nal/external approach, the transition from internal shock 
to external shock emissions takes place, typically, be- 
tween 50 and 100 seconds, precisely in the gap between 
BAT and XRT, at which time the flux levels of the in- 
ternal and external shock components are roughly equal. 
It may seem coincidental that the times for the tran- 
sitions from inter nal to external shocks are « (cf. 
IZhang et al.ll2006l ) . Swift observations of late-time X-ray 
flares interpreted as central engine activity mean, how- 
ever, that the internal shocks operate over a wide range 
of times, nullifying to some extent this coincidence. It 
may also appear coincidental that the flux level of the 
rising external-shock component is at about the same 
level as the earlier prompt internal shock component. 
But the large fraction (> 90%) of swept-up energy in 
the form of nonthermal hadrons, the factor 10 - 100 
more nonthermal energy in the shells before the radia- 
tive discharge compared to after, and the low efficiency 
~ 1 - 10% for transforming i nternal shock energy to 
hard X-rays and 7 rays J Belobor odovl[2000l : lKumajfl999l : 
iGuetta. Piran. fc Waxmani 120041. naturally leads to the 
expectation that energy radiated in afterglow phase is 
comparable to that of the prompt phase, so that the 
similar flux levels of the internal and external shock com- 
ponents may not be so surprising. 

We conclude that the proposed interpretation for the 
fast decays is in accord with the internal/external shock 
scenario for GRBs, provided that the external shock 
emissions make a significant contribution already in the 
early afterglow. 

5.3. External Shock Model 

A conceptually simpler though not generally accepted 
approach is to suppose that all GRB emissions are formed 
by a single impulsive event. In such a picture, the 7- 
ray pulses and X-ray fiares seen in GRB light curves 
are a consequence of interactio ns of the blast wave with 
a clumpy circumburst density (jPermer fc Mitmanlll999l . 



l2004f ). and there is a single blast wave that is energized 
only by the external forward/reverse shock pair. 

The central requirement for short timescale variability 
in this model is that the size scale of the inhom ogeneities 
(clouds) is < x/To (jPermer fc Mitman|[l999D . and that 
the blast wave remains thin in order to form a strong 
forward shock, which translates int o the condition tha t 
the blast wave spreading is < x/T^ (jPerme H [200l[2007l) . 
Even late-time X-ray fiares could be the result of a GRB 
blast wave intercepting clumps of mat erial if the blast- 
wave shell remains thin (|P erme 3 120071 ). A prediction of 
this model is that ranges of baryon loading, in p articular, 
mass- loaded dirty fireballs (jPermer et al.lll999l) . can ac- 
count for the differences between classical long duration 
GRBs and X-ray fiashes. 

This model h as been questioned on ground of efficiency 
(ISari fc Piran|[r9 97: Piran 2005), but the argument is in- 
correct (jPermer fc Mitmanii2004f) . The additional energy 
content in the prompt phase makes high efficiency for 
the external shock model, as for the internal shell model, 
less critical, improving the viability of both models. An- 
other concern is that late sharp flaring eve nts cannot be 
made by external shocks. The simulation of lZhang et al.l 
(120061) treats however a l arge cloud, and the treatment 
of iNakar fc GranotI (I2OO6I ) operates in the adiabatic self- 
similar phase, after which much of the blast-wave en- 
ergy has already been converted to radiation (this mis- 
take, and the failure to take into accou nt beaming factors 
for the interaction, is also made by lloka et all [2005). 
The brightest flares are made by portions of the blast 
wave that have not yet entered this phase, and which 
intercept a dense clump o f ma terial with a th ick column 
(jPermer fc Mitmanll 19991 ). See lPerrne 3 (f2007f) for a fuller 
discussion. 

Whether or not a "pure" external shock model for 
GRBs can be defended, the proposed model is already 
in accord with the internal/external scenario, as al- 
ready described . It is also interesting to note that 
iMcMahon et all (|2004D find that the afterglow 63 pa- 
rameter must be much weaker than in the prompt phase 
to explain light curves of certain smoothly varying GRB. 

5.4. Explanation for the Rapid Decay 

We have considered a physical mechanism to explain 
the rapid decays observed with Swift from GRBs in the 
late prompt/early afterglow phase that could operate 
in both internal/external and external shock scenarios. 
The number of physical models for this e ffect is not 
yet large; there is a preference recently (e.g . iNava et al 



'2007'; 'Panaitescu'2006': ' O'Brien et aI1l2006at IZhang et al 
2006; Liang et al. 2006) to quantify the effect and per- 
form empirical phcnomcnological studies rat her than to 
constr uct physical models. An exception is iPe'er etldl 
()2006l ). who consider scattering of GRB emission in a co- 
coon formed as the bl ast wave jet emerge s from the stellar 
photosphere (see also lPe'er et al.ir2006al ). The emergent 
photons can exhibit a steep temporal decay, similar to the 
GRB X-ray light curves observed with Swift XRT. This 
model applies strictly to a coUapsar scenario, where the 
shocks emerge from the stellar p hotosphere by forming 
a cocoon (e.g.. iRamirez-Ruiz et al. 2002), yet rapid de- 
cays also occur in short , hard GRBs (e.g., GRB 050724, 
IBarthelmv et al.ll2005ar ). which are thought to originate 
from coalescence events that lack stellar photospheres 
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and cocoons (for recent reviews of shor t hard GRBs, see 
lNakail[2?)07l : iLee fc R,amirez-R,nid[20?)1.3 

The curva t ure effect (iKuinar fc Panaitescul l200d : 
iDermeij 120041: iLiang et all l2006f) . which hmits the rate 
at which an impulsively illuminated shell can temporally 
decay, obviously plays an important role. The curvature 
limits apply to the case of an active central engine that 
suddenly shuts itself off, o r dramatically re duces its ac- 
tivity, as in the model of iProga fc Zhan"^ ([2006). For 
a spherical blast wave that is roughly uniform within 
the Doppler cone, the curvature limits would also ap- 
ply to internal or external shocks. The model presented 
here is in accord with the curvature limit, keeping in 
mind the uncertainty in se tting the zero of time (see also 
iLazzati fc Begelman|[200 6^. 

From a comparison of the idealized bolometric light 
curves in Figs. 5 -8 with the S wift BAT and XRT G RB 
light curves shown in, e.g., the lO'Brien et al.l (|2006D pa- 
per, it seems that the model offers adequate variety to 
explain the rapi d X-ray dec ays, and not simply an aver- 
age light curve ijO'Brien et al. 2006a') . 

5.5. Explanation for the Plateau Phase 

Focusing on the external shock component in the pro- 
posed scenario, about the time that Fermi processes can 
accelerate a large fraction of the particle energy to ultra- 
high energies, photohadronic losses and escape become 
important and cause a sudden loss of internal energy. 
The magnetic field energy, which is supported by the in- 
ternal energy, itself becomes smaller, thereby enhancing 
nonthermal particle escape and causing the internal en- 
ergy to collapse. At some point, the energy content in 
the blast wave reaches its minimum when only the orig- 
inal baryon load and a remnant magnetic field remains. 
Magnetic fiel d generation processes like the Weibel in- 
stability (e.g.. lNishikawa et al.ll200"6l ) could then amplify 
the blast wave magnetic field until plasma is swept-up 
and entrained, and the blast wave again commences to 
evolve adiabatically. 

The viability of this explanation for the plateau phase 
and the diversity of X-ray light curves discovered by Swift 
([Barthclmy et al. 2005) is suggested by the wide range 
of idealized light curve profiles calculated and plotted in 
Figs. 5-7, though more research will be required to 
quantify the underlying assumptions about the accelera- 
tion and particle radiation mechanisms that will allow fits 
to data. This explanation is quite different from the con- 
ventional explanation in the interna l/ external scenario. 
Ther e , delayed energy in jection (e.g.. lZhang fc MeszarosI 
120011: IZhang et al]|2006D . possibly signaled by the late- 
time X-ray flares, could harden the temporal decay. Yet 
plateau phases are often seen soon after a rapid decay, as 
in the cases cited in the Introduction, so the central en- 
gine activity has to be concealed while pow ering the light 
curve. An analysis of GRB 050713A by iGuetta et al.l 
(|2007t l starts with the conclusion that the behavior of 
GRB 050713A requires refreshed shocks. iToma et al.l 
(j2006') consider a patchy jet scenario where the delayed 
off-axis emissions become brighter than the on-axis emis- 
sions, but their model lacks a self-Compton component, 

^ This invites the speculation that the short hard class of GRBs 
are also sources of UHECRs, given appropriate parameters, in par- 
ticular, small variability times. 



and the afterglow X-ray component lasts much longer 
than found in detailed numerical c al culations with simi - 
lar parameters (jPermer et al.ll2000D . iFan fc Pira^ (|2006D 
try to model this behavior and assess the energy re- 
quirements. They claim that powering the hard after- 
glow phase does not make excessive demands on ener- 
getics, but do not display blast- wave solutions charac- 
teristic of the rapid decay and plateau phases discovered 
(|Nousek et al.ll2006l : IZhang et al.ll2006f l with Swift. 

5.6. Predictions 

Besides the steep decays and plateau phases in the X- 
ray light curves from GRBs, there are other direct ra- 
diative signatures for the hypothesis that UHECRs orig- 
inate from GRBs. In all cases where photopion losses are 
important, the blast wave is assumed to be radiating in 
the fast-cooling limit, so that the I'Fi, indices a = 1/2 
and 6=1— (s/2), where s is the electron injection in- 
dex (with s w 2.2 for shock Fermi models). If we accept 
that accelerated electrons radiatively cool without be- 
ing subsequently accelerated (contrary, incidentally, to 
the point made earlie r with respect to the concerns of 
iGhisellini et al.l ()2000l )). then already a careful spectral 
analysis for the variation of Cpk predicted by external 
shock emissions as it sweeps through the Swift XRT de- 
tector could see if GRB blast waves evolve within a few 
hundred seconds to a condition where the eg and pa- 
rameters are each > 10%. 

A prediction of this model is hadronic 7-ray light con- 
sisting of proton synchrotron, photopion and secondary 
photohadronic decay radiations that cascade to low 
enough energies that internal and diffuse cosmic 77 at- 
tenuation processes become negligible. This component 
varies differently from the X-ray lepton synchrotron com- 
ponent and can appear much after the prompt emissions. 
The hadronic emission component could then possibly 
account for the del ayed 7-ray emission in GRB 940217 
(iHurlev et al.l [l99l and GRB 941017 (iGonzalez et all 
l2003f ). likelv in addition to SSC emission. The inde- 
pendently varying hadronic emission component is a 
prediction for GLAST , though a leptonic model (e.g., 
iPe'er fc WaxmanI I2004D for delayed anomalous compo- 
nents that includes Compton-scattered forward- and 
reverse-shock emissions, could also make a delayed com- 
ponent as observed from GRB 941017. 

Ultra-high energy neutrinos are made by the release 
of blast wave energy though this mechanism, and we 
predict the detection of a delayed neutrino flux at high 
energies just about the time of the rapid decays. Un- 
fortunately, the > 10"'^^ eV neutrinos from the photo- 
hadronic processes in GRB blast waves are at a difficult 
energy for IceCube to detect, but the > 10^^ eV neutri- 
nos could be detected in sideways showers with Auger, 
or in balloon-borne detectors such as ANITA. Besides 
the UHE neutrino component, PeV neutrinos detectable 
with IceCube will be made during the prompt phase in 
the internal shcU/coUapsar model in large fluence GRBs 
when the shocked shell Lorentz factors ~ 100 - 200, 
while bei ng optically thick to 77 attenuation at GeV 
ener gies JPermer fc AtovanI [2003: iRazzaaue et al.l 12004 
Guetta^ei al. 2004). 

For parameters that make a highly radiative blast 
wave phase, ep^syn for proton synchrotron emissions from 
UHECRs in GRBs wfll be found in the TeV range. To 
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avoid attenuation by the diffuse intergalactic infrared ra- 
diation fields, this emission could only be detected from 
low redshift, z < 0.4, GRBs, which can be measured with 
low-threshold air Cherenkov telescopes such as MAGIC. 
The photon flux at these energies is very low, so only a 
few > 100 GeV photons could be detected with a direct 
photon detection telescope such as GLAST. 

6. SUMMARY AND CONCLUSIONS 

By carrying through an analysis of a complete lep- 
tonic/hadronic GRB blast- wave model, a region of pa- 
rameter space was found for GRB external shocks where, 
even in the simplest constant density blastwave model, 
photohadronic processes have important effects on blast- 
wave evolution. This occurs for the observer a few hun- 
dred seconds after the start of the GRB for circumburst 
medium densities no > 10^ cm'^^ lOO < To < 300, 
£■54 ~ 1 - 10, and ee,eB ^ 10%. For these parameters, 
the time scales for photopion losses and direct escape 
become comparable to the available time in the early af- 
terglow, with observable consequences. 

The GRB blast wave is predicted to evolve toward the 
strong cooling regime in the early afterglow, leading to 
a strong discharge of energy through photohadronic pro- 
cesses and direct escape of UHECR ions, causing the de- 
clines in the X-ray light curves observed by Swift. The 
change of Cpfc with time may indicate the evolution of the 
radiative regime during the prompt and early afterglow 
phases, though this component can be concealed in an in- 
ternal/external shock model by internal shock emissions. 



Hadronic radiations consisting of high-energy neutrinos, 
cascade 7 rays, and escaping cosmic-ray neutrons could 
be important to > 10"* s after the GRB, producing a 
component that varies, in general, differently from the 
X-ray lepton sychrotron emission. As a prediction for 
GLAST, this component should be bright at GeV ener- 
gies just before the times that the rapid X-ray declines 
are observed. A related flux of > 10^^ eV neutrinos is 
made at this time. 

The release of energy in the form of neutrons or escap- 
ing ions in the early afterglow phase means that there 
was much more energy in the GRB blast wave during the 
prompt phase. Thus the efficiency of 7-ray production 
during the prompt phase is low, which ameliorates effi- 
ciency concerns in both external or internal shock mod- 
els. Hadronic processes imprint GRB light curves and 
the UHECR spectrum with distinctive signatures. The 
rapid X-ray declines and plateau phases discovered with 
Swift are argued here to be just such signatures. 
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